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Low-mass star formation in R Coronae Australis: Observations of 
organic molecules with the APEX telescope 
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Abstract. This paper presents new APEX submillimetre molecular line observations of three low-mass protostars, IRS7A, 
IRS7B, and IRAS32, in the R Coronae Australis molecular cloud complex. The molecular excitation analysis is performed 
using a statistical equilibrium radiative transfer code. The derived beam averaged fractional abundances vary by less than 
a factor of two among the three sources, except those of IFCO and CH3OH, which show differences of about an order of 
magnitude. The molecular abundances are similar to those typically found in other star-forming regions in the Galaxy, such as 
the p Oph and Perseus molecular clouds. There is a marked difference in the kinetic temperatures derived for the protobinary 
source IRS7 from HiCO (~ 40-60 K) and CH3OH («20K), possibly indicating a difference in origin of the emission from 
these two molecules. 
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1. Introduction 

In the earliest, deeply-embedded stage, low-mass protostars 
are surrounded by a collapsing envelope and a circumstellar 
disk. In this phase, the initial chemical composition of the disk 
is determined. The chemistry of molecules in the circumstel- 
lar envelope is likely to be directly reflected in the molecular 
composition in the disk, which forms the building blocks from 
which planets may form. An interesting aspect is the search 
for low-mass 'hot cores', i.e., regions where grain ice mantles 
are liberated in the envel ope close to the centra l protostar be- 
cause of its heating (e.g. , ICeccarelli et alJl200fi ISchoier et ail 
2002i lCazaux et alJ2003l) . The search for hot cores around low- 
mass stars is not yet conclusive, mainly due to limitations in 
the available observational data. Also, outflows are likely to 
play an important role by affecting many of th e same species 
thoug ht to be sign posts of hot cores (e.g., I.Tc>rgensen et alJ 
2005). Moreover, the chemistry of traditional hot-core tracers 
such as H?CO may be affected by the freeze-out of CO (e.g., 
ISchoier et all2004l) . which is why combined studies of chemi- 
cally related molecular species are important. 

The R Coronae A ustralis dark cloud is a nearby (170pc, 
Knude & H0g 1998) region of active star formation (e.g., 
BrowtJ ll987t IChini et aljkoolt iNutter et alJl2005fc . Molecular 



have attained little attention, mainly due to its location on the 
southern hemisphere. Previous studies using single-dish tele- 
scopes with limited spatial resolution, which are mainly sensi- 
tive to the chemistry of the outer cold regions of the envelopes, 
have been biased towards lower excitation lines. A wealth 
of molecular rotational transitions lying in the 279 - 381 GHz 
window are accessible with APEX 1 , and can be used to con- 
strain the physics, and in particular, the chemistry in the 
dense (> 10 6 cm 3 ) and warm (> 30 K) material in the envelopes. 
Observations of protostellar objects in Corona Australis are 
also interesting for testing the importance of the surroundings 
on the formation of protostars, and will complem ent previ- 
ous extensive stud ies in Perseus and p Oph Jj0rgensen et alJ 
120021 l2004i 120051) . It is not clear how the properties of the 
pare ntal cloud affect the prop erties of the emerging protostar 
(seellavaward hana et aljEoOll) . Likewise, the properties of the 
parental cloud may reflect the chemical evolution of the pro- 
tostar, for example, in the degree of freeze-out and the 'initial 
abundances ' . 

The first step in a project aimed at obtaining an exten- 
sive data base of molecular line emission towards a large sam- 
ple of protostars in the southern hemisphere, with a focus on 
the Corona Australis and Chameleon molecular clouds, is pre- 



line observations of protostellar objects in R Coronae Australis 
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1 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory 
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Fig. 1. New observations of molecular line emission using the APEX telescope. The velocity resolution is 0.2 km s 1 in all the 
spectra. 



sented in this Letter. It is crucial to utilize APEX to build com- 
parable samples of well-studied low-mass protostars for future 
ALMA observations 

2. Observations 

Observations of sub-millimetre molecular line emission to- 
wards the R Coronae Australis molecular cloud complex were 
performed September 4-6, 2005, using the 12 m APEX tele- 
scope equipped with an SIS receiver (APEX-2a). Typical sys- 
tem temperatures were in the range 140- 380 K (DSB). The 
1 GHz FFTS with 8192 spectral channels was used as a back- 
end, providing a spectral resolution of 122kHz (0.1 kms~'). 
The observations were carried out using a position-switching 
mode, with the reference position located at +2' in azimuth. 



The data were reduced in a standard way by removing, in most 
cases, a low order (< 3) polynomial baseline and the raw spec- 
tra, stored in the T* K scale, converted to main-beam brightness 
temperature using T m \, = T^/rj m \,. T* K is the antenna tempera- 
ture corrected for atmospheric attenuation using the chopper- 
wheel method, and r]^ is the main-beam efficiency. A value of 
?7mb = 0.7 was assumed for all spectral settings. Regular point- 
ing checks were made on strong CO sources and typically 
found to be consistent within « 3". The FWHM of the main 
beam is 18"at 345 GHz. The uncertainty (including errors in 
pointing, calibration, and t]^) in the absolute intensity scale 
(main-beam brightness) is estimated to be about +15%. 

The observed spectra for the protostellar objects 
IRS7A (a 20 oo = 19:01:55.3, $2000 = -36:57:21), IRS7B 
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Fig. 2. New observations of CH3OH line emission using the 
APEX telescope. The velocity resolution is 0.5kms _1 in both 
spectra. Note that the 7o — > 60 E + line is blended with the 
CH3OH 4q — > 3 1 E~ line from the upper sideband. 



(0:2000 = 19:01:56.4, £ 20 oo = -36:57:27), and IRAS32 
feooo = 19:02:58.7, #2000 = -37:07:34) are presented in 
Figs. H an d ID (the positions are obtained from IRAC4 8 /im 
pipeline data in the Spitzer Space Telescope archive: AOR ID 
3650816). IRS7B is identified as a cl ass protostar, wh ereas 
IRS7A and IRAS32 are of class I (see lNutter etai]l2005l and 
references therein). 

Results of Gaussian fits to the spectra are reported in 
Table ^ The systemic (LSR) velocities obtained from the 
Gaussian fits are 5.3, 5.6, and 5.9 km s" 1 (LSR) for IRS7A, 
IRS7B, and IRAS32, respectively. The spectra are generally 
well represented by a single Gaussian, with the major excep- 
tion being the observed HCN J — 4 — > 3 line emission towards 
IRS7A and IRS7B. In these two cases, not only is the HCN line 
significantly broader, compared to other molecular lines, but 
also a central absorption dip is present, at least in the IRS7A 
spectrum. The HCN spectra presente d here are very s i milar t o 
those of HCO + 7 = 4^3 obtained by And erson etail dl997ah . 
They interpret the emission as arising in a large rotating (cir- 
cumbinary) disk, and the self-absorption feature is caused by 
absorption in a dense and quiescent foreground of the cloud. In 
addition, for the two IRS7 sources, our observed C n O and, in 
particular, C ls O, J — 3 — > 2 spectra have weak wings, possibly 
related to the known molecular outflow in this region, and are 
somewhat better reproduced by a sum of two Gaussians. 

3. Excitation analysis 

Presently, there is a lack of observational data required to carry 
out a detailed modelling of t he dust emission, m ainly due to 
limited spatial resolution (see N utter et alJl2005l) . This is usu- 
ally the preferred method to determine the physical structure, 
such as density and temperature, of material present around 
protostars that is required for a detailed analysis of its molecu- 
lar line emission ("e.g.. l.TOrgensen et all20 02). 

A first estimate of the excitation conditions can be obtained 
assuming the lines to be emanating from the same homoge- 
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a /= f r mb rfv is the integrated intensity in Kkms~'. T mb is the peak 
intensity (in K) and Av the FWHM (in kms~') from Gaussian fits. A 
colon (:) marks an uncertain value due to low a S/N-ratio or a signifi- 
cant departure from a Gaussian line profile. 

b Blended with the CH 3 OH 4 — * 3] E~ line from the upper sideband. 
The integrated intensity refers to the sum of both lines. 



neous region. We have adopted RADEX, a statistical equi- 
librium radiative transfer code that uses an escape probabil- 
ity formalism to calculate the excitation and emission from 
the protostellar cores in our sample. RADEX is comparable 
to the LVG method and provides a useful tool for rapidly 
analysing a large set of observational data providing constraints 
on physical conditions, such as density and kinetic tempera- 
ture (Jansen et al. 1994; van der Tak et al., in prep.). An on- 
line version of RADEX, together with relevant molecular data 
used in the excitation analysis and summarized in lSchoier et alJ 
ll2005h . is made publicly available through the Leiden Atomic 
and Molecular Database (LAMDA) 2 . 

In the excitation analysis, the number density of H2 
molecules was fixed to 3 x 10 6 cirT 3 , typical of the warmer, 
more dense regions of the condensed cores. From the two 
lines of H2CO, the kinetic temperatures were estimated to be 
*60K and ^40K for IRS7A and IRS7B, respectively. The 
corresponding column densities derived were 2xl0 13 cm~ 2 
(IRS7A) and 3.5 x 10 13 cirT 2 (IRS7B). In the case of CH 3 OH, 
a full analysis using a ^-statistic was used, and the kinetic 
temperatures derived were 23+3 K and 22+3 K for IRS7A and 
IRS7B, respectively. The column densities of CH3OH derived 
were 1.4±0.4x 10 14 cm 2 (IRS7A) and 1.8+0.4 x 10 14 cirT 2 
(IRS7B). The substantial difference in the derived kinetic tem- 
peratures from H2CO and CH3OH could possibly indicate dif- 
ferent origins for these traditional hot-core tracers. 



2 http : //www. strw. leidenuniv.nl/~moldata 
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Table 2. Model results. 
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A colon (:) marks an adopted value used to set the density scale. 



we get for H2CO (~ 40-60 K), suggests the impact of outflows 
in regulating the chemistry around IRS7. These t emperatures 
are s i milar to those obtained for NGC-IRAS4 JMaret et alJ 
l2004 12005b, but lower than those for IRAS 16293-2422 
Jvan Dishoec k et a"fll 19951) . However, in contrast to the NGC- 
IRAS4A binary in particular, IRS7 does not show CH 3 OH 
lines that are broader than other observed molecular lines. 
Nonetheless, a large scale outflow thought to be associated 
with IRS7 has been observed in several other molecular species 
jLevreaultll988UAnderson et alJl997bh . 

To make further progress, in particular to solidify the claims 
regarding the chemistry in the direct neighbourhood of the ob- 
served protostars, more detailed knowledge of the underlying 
physical structure is required. Such work is in progress and 
involves mid-IR observations from large ground-based tele- 
scopes and the Spitzer space telescope, and interferometric 
mm/sub-mm observations from ATCA and the SMA. 



To determine the column densities for the remaining molec- 
ular species, we adopted a kinetic temperature of 40 K. The 
column densities obtained for C 34 S in the case of IRS7A 
(lxl0 13 crrT 2 ) and IRS7B (1.5 x 10 13 cm" 2 ) were used to- 
gethe r with a 'standard' fr actional abundance (5xl0~ n ) of 
C 34 S ( J0rgens en et alJl2004l see Table [5J to estimate the total 
H2 column densities. In the case of IRS7A, we obtain a H2 col- 
umn density of 2 x 10 23 cm" 2 , and for IRS7B, a slightly higher 
value of 3 x 10 23 cm" 2 . For IRAS32, no C 34 S 7 = 7^6 line 
emission was detected, and instead the main isotope (with an 
adopted fractional abundance of 1 x 10 -9 ) was used to estimate 
the H2 column density to 3 x 10 22 cm" 2 . The derived fractional 
abundances for all molecular species observed are reported in 
Table|2 

4. Results and discussion 

The excitation analysis gives results (see Table|2j for the abun- 
dances of common molecules such as CO, CS, and HNC, and 
for more complex species, such as H2CO and CH3OH, that are 
consistent with results fro m other star-forming regions, such 
as p Oph and Perseus ttorgensen et all2002ll2004ll2005l> . The 
derived fractional abundances vary by less than a factor of two 
among the three sources for the majority of molecular species. 
In the present analysis, it is not possible to correlate the abun- 
dances with the identified protostellar class. It is interesting to 
note that the abundances of both H2CO and CH3OH vary by 
about an order of magnitude among the three sample sources, 
possibly indicating a difference in central heating source or im- 
pact of outflows. We find the abundances of these molecules to 
be particularly low in IRAS32. 

IRS7A and IRS7B form a protobinary system with a sepa- 
ration of about 18" (3000 AU). As such, RCrA-IRS7 is also 
a cousin of the well-studied IRAS 16293-2422 and NGC- 
IRAS4 b inaries in term s of the total luminosity. The fact that 
Jeirgense n et alJ (12005) only detected CH3OH emission from 
known outflow sources in their large survey of protostellar 
cores, together with the low kinetic temperature (^20K) we 
derive from our observed CH3OH emission in contrast to what 
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